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ABSTRACT 

We have detected sporadic, bright, short-duration radio pulses from PSR J0901—4624. These pulses 
are emitted simultaneously with persistent, periodic emission that dominates the flux density when 
averaging over many periods of the pulsar. The bright pulses have energies that are consistent with 
a power-law distribution. The integrated profile of PSR J0901—4624 is highly polarized and shows 
four distinct components. The bright pulses appear to originate near the magnetic pole of the pulsar 
and have polarization properties unlike that of the underlying emission at the same pulse phase. We 
conclude that the bright pulses represent a secondary giant-micropulse emission process, possibly from 
a different region in the pulsar magnetosphere. 

Subject headings: pulsars: general — pulsarsiindividual (PSR J0901—4624) —stars:neutron 


1. INTRODUCTION 

Since the discovery of pulsars, the phenomenology as¬ 
sociated with their radio emission has been a fertile 
source for studying coherent-emission processes in mag¬ 
netized plasma. It is clear that although the time- 
integrated profiles of radio pulsars are constant or nearly 
constant, individual pulses that go into creating the inte¬ 
grated profile are highly variable both in total intensity 
and in their polarization. 

Historically, there has been an understanding that the 
fluence of individual puls es follows either a log-normal 
or Gaussian distribution (iCa irn s. Johnston &: Dasll200lL 
120031 : iBurke-Spolaor et alJ 120121 : [Shannon et alTl2014H 
with relatively few pulses significantly stronger than 
the mean. We will refer to this distribution as the 
normal-mode of emission. The main exception to this 
has been the giant pulses o bserve from the Crab pulsar 
(|Staelin fc ReifensteinlllM^ . which are extremely bright 
(> 1 MJy) and na rrow (<C 10 ns) burs ts emitted at cer¬ 
tain pulse phases (jHankins et al.l[200ll . 

More recently, a sub-set of pulsars have been iden¬ 
tified that show highly erratic behaviour in their sin¬ 
gle pul ses. The giant-pulse mo de was seen in other 
young (iJohnston fc Romai^ 200311 and millisecond pul¬ 
sars (iCognard et al.l Il99a iRomani fc JohnstonI 120011 : 
iKnight et al.l 1200611 . Several young, energetic pulsars 
show single-pulse variability similar to the giant-pulse 
mode, except that the bursts have us rather than n s 
widths (IJohnston et ahll^OOlUJohnston fc Romanill^002ll . 
which has been referred to as the giant-micropulse emis¬ 
sion. 

Emission from less energetic pulsars and millisecond 
pulsars also shows phase dependent variability, including 
power-law statistics on the leading or trailing edges. 

One of the best-studied slow pulsars, PSR B1133-I-16 
shows, at high frequencies > 4 GHz) in dividual pulses 
with power-law statistics, of ^ 2 ms width ((Kramer et al.l 
120031 : iKrzeszowski et all 1201411 . The bright pulses orig- 
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inate predominantly from the trailing edge of the first 
component of the double-cone profile. At low frequen- 
cies, the bright pul ses originate on both components 
((Kramer et al.l 1200311 . The brightest pulses from the 
millisecond pulsar PSR J1713-I-0747 originate from the 
traili ng edge of the main component ((Shannon fc CordesI 
l2012fl . However these pulses show a lognormal energy 
distribution. In contrast, the brightest pulses in PSR 
J0437—4715 originate from the center of the profile. 
These stronger pulses h ave fractionally highe r polariza¬ 
tion than weaker pulses ((Oslowski et al.ll2014ll. _ 

Ro tating Radio Transients (RRATs. lMcLaughlin et al.l 
120061 1 emit single pulses extremely sporadically, with 
pulses emitted every few minutes or longer. While 
some of the pulsars in this group also ha ve un¬ 
derly ing normal-mode emission, some do not ((Ridlevl 
l2010fl . Further complicating the picture, several 
pulsars appear to have a co mbination of normal¬ 
mode, and RRAT-like behaviour (IWelteyrede et al.ll2006t 
IWeltevrede. Johnston fc EsDinozall2011ll . 

Here we report the discovery of a new example of spo¬ 
radic, bright emission from PSR J0901—4624. Bright 
pulses from J0901—4624 show a power-law distribution, 
which are weaker and broader than giant pulses but occur 
more frequently than giant micropulses identified from 
previous pulsars. In addition, we compare the bright 
pulses of J0901—4624 to other modes of pulsar emission. 

2. OBSERVATIONS 

During a monitoring programme of ^ 180 young pul¬ 
sars at the 64-meter Parkes radio tel escope in support 
of th e Fermi Space Telescope mission (IWeltevrede et al.l 
l2010f). using a rea l-time transient search algorithm 
((Barsdell et al.l (2012h . we identified occasional (once per 
20 s) bright (S/N > 20), narrow (< 500 fis) pulses be¬ 
ing emitted from PSR J0901—4624. With a rotation pe¬ 
riod of 442 ms, a dispersion measure of 199 pccm“^, 
and an inferred spin-down energy of 4.0 x 10^"'^ ergs“^, 
PSR J0901—4624 has median properties for the sample 
selected for high spin-down energy. 
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To examine the bright pulses in more detail, we ob¬ 
served PSR J0901—4624 again with the Parkes telescope 
on six separate occasions, of 30 — 60 min in duration, as 
listed in Table [H All observations were made using the 
21-cm multibeam system at a central frequency close to 
1400 MHz. The system-equivalent flux density of the 
system is ^ 35 Jy on cold regions of the sky. 

Data were recorded using three backend systems oper¬ 
ating in parallel. The first backend (PDFB4) produced 
a period-averaged digital-filterbank output of 256 MHz 
bandwidth, centered at 1369 MHz, comprising 1024 fre¬ 
quency channels for each of 1024 phase bins across the 
pulse period for both the two auto correlations and 
the real and imaginary parts of the cross correlations 
of the feed probes. Data were folded at the topocen- 
tric pulse period and accumulated for 30 s and recorded 
to disk. The second system (PDFB3) produced fast- 
sampled digital-filterbank output with 256 MHz band¬ 
width centered at 1369 MHz, with 512 frequency chan¬ 
nels per polarization sampled every 256 ^s. Finally, on 
May 25 we used a baseband recording system (CASPSR), 
with 400 MHz bandwidth centered at 1382 MHz in or¬ 
der to resolve short-duration signals and mitigate intra¬ 
channel dispersion smearing in observations with the 

filterban k-spectrometer system s. _ The DSPSR software 

package l|van Straten fc BailesI 1201111 was used to sub¬ 
divide the data into single rotations of the pulsar for the 
PDFB3 and CASPSR datasets. In the case of CASPSR, 
DSPSR was also used to coherently disperse the data. 

A calibration signal, injected into the feed at an angle 
of 45° with respect to the probes, was recorded before 
each observation. This signal was used to determine the 
relative gain and phase of the two polarization probes. 

Data editing and ca libration were conducted using the 
PSRC HIVE package (jHotan. van Straten fc Mancheste'^ 
l2004fl . In brief, after removal of radio-frequency interfer¬ 
ence in both the frequency and time domains, the data 
were gain and polarization calibrated and summed in fre¬ 
quency and time to produce a pulse profile for each pulsar 
observed. Data were flux calibrated using observations 
of the radio galaxy Hydra A m ade as part of the Parke s 
Pulsar Timing Array project (jManchester et aH 1201311 . 
Details of the observations can be found in Table [TJ 

3. INTEGRATED PROEILE 

In the top panel of Fig [l] we show the full- 
Stokes profile and linear polarization position angle for 
PSR J0901—4624, formed from co-adding all of our ob¬ 
servations. The profile consists of four distinct compo¬ 
nents with a total width of 27°. The geometric cen¬ 
ter of the profile occurs between the two inner compo¬ 
nents, which indicates the pulsar has double-conal struc¬ 
ture. The pulsar is highly linearly polarized, which is 
typical for young, high-A pulsars. PSR J090I—4624 
also shows a high degree of circular polarizat i on in 
the trailing component (iJohnston fc Weisberd l2006t 
IWeltevrede fc Johnstonlf2008fl : it is common for the trail¬ 
ing component to have the do minant amplitude and 
signif icant circular polarization (|.Iohnston fc Weisberd 

Ho^. 

The position-angle swing is well described by the 
rotating-vector model with the best-fitting model giving 
an inflexion point located 3.9° after the geometric cen- 
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Fig. 1.— Upper panel: Integrated profile of 
PSR J0901—4624 in full Stokes at 1369 MHz The lower 
plot shows the total intensity (black line), linear polar¬ 
ization (red line) and circular polarization (blue line) as a 
function of pulse longitude. The top panel shows the po¬ 
sition angle of the linearly polarized radiation and runs 
between —90° and -1-90°. The green line is the best- 
ftting model of the beam geometry to the position-angle 
data. Phase zero is chosen to be the geometric center of 
the profile. The dominant features in the profile are four 
components centered at longitudes of —12°, —6°, 4°, and 
9° relative to the geometric center. Lower panel: Profile 
formed from only the bright (S/N > 5) pulses. Bright 
pulses originate primarily from the central component 
and have a width of 1-2 bins of phase resolution. The 
RVM model is the same as in the upper panel. 
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TABLE 1 

Observation Details 


Date 

(MJD) 

N 

N5 

s. 

(mJy) 

May 17 

8148 

32 

0.66 

May 25 

8156 

33 

0.69 

July 25 

7034 

72 

0.95 

July 27 

2197 

17 

0.85 

Aug 19 

4085 

13 

0.63 

Sep 28 

4085 

5 

0.51 


Note. — For each epoch, we list 
the N, number of rotations of the 
pulsar observed; N 5 , the number of 
5 — (7 pulses detected, and Su, the 
mean observed flux density. 

ter of the profile. This can be used to infer an emission 
height of ^700 km, or about 3% of the light-cylinder ra¬ 
dius. Using the method of Rookyard et al. (2014) we 
determine a best-fitting value for the inclination of the 
rotation axis to the magnetic axis of a ~50° and an im¬ 
pact parameter of /3 ~2°. 

Observations of the pulsar have also been made at 
3.1 GHz and 0.7 GHz. There is little evolution of the 
profile over this frequency range. 

4. INDIVIDUAL PULSES 

In total, over the six observations we observed ^34,000 
rotations of the pulsar. Approximately ^400 pulses with 
a S/N > 5 were identified. All of these bright pulses 
had narrow width ( < 1 ms). The daily variation in the 
rate of bright pulses, as shown in Table [U is caused by 
scintillation and the distribution of pulse energies (see 
section El). After accounting for these effects, consistent 
numbers of bright pulses were detected on each day. 

Based on observations made with the GASPSR coher¬ 
ent dedispersion system, the full-width half maximum 
of the bright pulses varies between 130 /is and 400 ^s, 
which is less than 0.1% of the pulse period. Two of the 
brightest pulses from these observations are displayed in 
Figure E] The location of the bright pulses varies in lo¬ 
cation by approximately 4° in pulse phase (1% of the 
period). Occasionally bright pulses are seen from the 
two components that flank the central component. The 
time between bright pulses is consistent with a Poisson 
random process. 

In the lower panel of Figure [H we show the average 
profile formed by combining only the brightest pulses. 
These pulses greatly enhance the emission at the center 
of the integrated profile. We note that the rarity of bright 
pulses results in very little overall effect on the integrated 
profile. In addition, we see the bright pulses are emit¬ 
ted simultaneously with the underlying, fainter emission. 
The normal emission observed during the bright pulses 
(‘S'peak. bright ^ 50 mjy) is modestly brighter during all 
of the pulses (S'peak.aii ^ 35 mJy) suggesting that there 
may be a connection between the two modes of emission. 

4.1. Polarization Properties 

The integrated profile of PSR J0901—4624 is highly 
linearly polarized (top panel of Figure 1). This im¬ 
plies that one orthogonal mode dominates the emission 
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Fig. 2.— Two bright pulses taken from our coherently 
dedispersed observations with GASPSR. The data have 
been binned to provide 32 p,s resolution. Panel a Full 
rotation of PSR J0901—4624. Panels b zoom-in of pulse 
displayed in Panel a. Panel c: Second bright pulse. 

((Johnston et al.l 120051) and that the single pulses should 
also be highly polarized. For the normal emission this 
is indeed the case. We occasionally detect bright pulses 
from component 2 and 4 which are highly polarized (see 
lower panel of Figure 1). 

In contrast, the bright emission has very different po¬ 
larization as can be seen in Figure I. In the centre of 
the profile, where the bright emission originates, almost 
no polarization is seen and what little polarization there 
is has a different polarization-angle swing to that of the 
normal emission. This is unusual behaviour and hints at 
a different origin for the normal and bright emissions. 


5. FLUX DISTRIBUTION 

We analyzed the flux distribution using the single-pulse 
PDFB4 data. We identified bright pulses by first iden¬ 
tifying the brightest bin in a window of 14° width in 
pulse phase centered on the third component of the av¬ 
erage profile. We then calculated three signal-to-noise 
(S/N) ratios using an estimate of the off pulse noise: the 
first from only the brightest bin, and the second and 
third from the sum of the brightest and its two adjacent 
bins. We only searched for narrow pulses after a first in¬ 
spection of the data revealed that the bright pulses were 
either unresolved (< 1 bin of pulse phase for the lower 
time resolution PDFB4 data) or marginally resolved. Of 
these three, we report the integrated pulse energy with 
the largest S/N. To assess the expected distribution of 
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noise, we repeated the same algorithm using a control 
sample of the same width in a region of pulse phase in 
which no pulsar emission occurs. In Figure |3l we show 
that the histograms of pulse energy in both the bright 
pulse window and the off-pulse window. 

Because of the relatively low signal-to-noise ratio, it 
is necessary to deconvolve the noise distribution when 
modelling the pulse -energy distribution. We used tech¬ 
niques described in [Shannon et al.l (1201411 to deconvolve 
the noise and model the pulse-energy distribution. The 
measured pulse energy is the convolution of the noise 
Pn and the intrinsic energy pj: 

Pe{E)= j dE'pN{E')pi{E-E'). (1) 

The noise distribution pjq is non-trivial to calculate 
because the algorithm used to identify significant pulses 
results in a non-Gaussian statistics, and refractive scin¬ 
tillation modulates the S/N at each epoch. For a single 
observing epoch, in which we assume that the scintilla¬ 
tion does not significantly affect the pulse energy, by the 
central limit theorem, the off-pulse samples are expected 
to be Gaussian distributed. The probability of getting at 
least one measurement greater than some value E in N 
trial samples is the converse of the probability of having 
all values less than E 


P{E > E') 


1 

2 


1 -I- erf 



( 2 ) 



Pulse Energy (mJy ms) 

Fig. 3. — Pulse-energy distribution for 
PSR J0901—4624. The thin solid line labelled N 
shows a model for the noise. The thin-dashed line 
labelled S shows the best-htting power-law model for 
the pulse energy distribution. The thick solid line 
(labelled S + N) shows the convolved pulse-energy 
distribution. 


The probability density is the derivative of Equation 


N 

. ^ 1 


N-1 

x‘^ 

y 

1 -I- erf 

VV^CT/. 

exp 



Because we have averaged together multiple days, af¬ 
ter correcting for pulse flux variations caused by scintilla¬ 
tion, the noise distribution deviates from theoretical ex¬ 
pectation. We therefore fitted an analytic function to the 
off-pulse noise distribution, shown as the thin, solid line 
in Figure 131 with N and a allowed to be free parameters, 
and we then deconvolved the pulse-energy distribution 
using this function. 

For the pulse energy distribution we c onsidered both 
log-no rmal distribution (see Equation 7 of iShannon et al.l 
(|2014f) l and power-law distribution 

pi{E) = AE-^ eM-Ec/E) {E < E,) 

= AE-^ {E > E,), (4) 

where /3 is the spectral index, Ec is a low-energy expo¬ 
nential cutoff, and A is a normalising factor. 

We find that models with power-law distributions pro¬ 
duced better agreement than models with log-normal dis¬ 
tributions. In particular the best-fitting log-normal mod¬ 
els do not reproduce the observed high-energy tail. For 
the best-fitting power-law model, the energy cut-off is 
Ec ~ 0.1 mJy ms and the spectral index is a = 2.28(4). 
The spectral index is comparable to those measured in 
other power-law emitting pulsars (iKnight et al.ll200^ . 

6. DISGUSSION 

One of the major open questions in pulsar emission is 
how timescales are set for pulse emission. Giant pulses 


have timescale of n s, whereas the intermittent pulsars 
(|Kramer et al.l 1^0611 can remain in quiescent states for 
years between active states. In between these two ex¬ 
tremes are microstructure, single pulse drifting, nulling, 
mode-changing and RRAT emission. There is clearly a 
continuum of temporal regimes, but it is important to 
assess how these variations relate to other pulsar prop¬ 
erties such ages, viewing and beaming geometries, and 
magnetic field strengths. PSR J0901—4624 shows a sec¬ 
ondary emission mode in which bright narrow pulses are 
frequently emitted from the central component. 

Another important differentiator between emission 
processes is the distribution of the energies of individual 
pulses (see ICairns. Johnston fc Das|[200lL and refer¬ 
ences within). Normal-mode emission (from slow or mil¬ 
lisecond pulsars) is observed t o have either log-normal or 
gaussian energy distr ibutions (|Burke-SDolaor et al.ll2012l : 
iShannon et all [20141 1. The situation with RRAT-mode 
emission is much less clear. The best-studied object 
with RRAT emission, PSR J1819—1458 shows sporadic 
pulses that are very sirn ilar to normal-mode emission 
(iKarastergiou et all 1200911 . Howe ver, other RR ATs ap¬ 
pear to have power-law emission (jRidlevir20ir)H . In con¬ 
trast, giant pulses, giant micropulses of Vela, the trail¬ 
ing bright pulses in PSR B1706—44, and the emission of 
PSR J0901—4624 follow power-law distributions. 

The width of the bright pulses in PSR J0901—4624 
seem to lie between the giant pulses (with widths less 
than 10“®P) and those of the RRATs (with widths 
greater than 10“^P). The jitter in the arrival times is 
also greater than in total pulse phase giant-pulse emis¬ 
sion, less than RRAT-mode emission, and similar to gi¬ 
ant micropulse emission. The widths of the pulses can be 
connected to the physical scales of the emission region, 
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with shorter duration pulses implying physically smaller 
emission regions. 

The location of the bright emission also can be used to 
distinguish emission modes. Giant-pulse emission is as¬ 
sociated with gamma-ray emission and like ly arises from 
the outer magnetosphere (lAbdo et al.ll20ir)fl . In contrast, 
the bright pulses observed from Vela, PSRs B1706—44, 
B1046—58, B1133-I-16, and J1713-I-0747 are not associ¬ 
ated with higher energy emission but arise at the edges 
of the normal-mode emission. Like PSR J0901—4624, 
the bright pulses from PSR J0437—4715 arise from close 
to the center of the profile; however these pulses show 
higher polarization than weaker pulses. 

In at least some of the RRATs, emission appears to 
come from conventional core- cone emission morphology 
across the whole polar cap (IKarastergiou et all 120091 
IKeane et al.ll201lll . In PSR J0901—4624, the bright pul¬ 
sars seem to originate close to the center of the profile 
and hence the magnetic pole. There is no detected high- 
energy emission for PSR J0901—4624, so we do not know 
if its bright pulses are coincident with a high-energy emit¬ 
ting region. 

Finally, the polarization properties of these emission 
modes are different. Both giant-pulse emission in the 
Crab (jHankins et al.l 1200311 and giant micropulse emis¬ 


sion i n the Vela pulsar ( Kramer. Johnston fc van StratenI 
I 2 OOI are highly polarised. For the RRAT 
PSR J1819—1458, the polarization properties are 
similar to that seen in single pulse emission; the degree 
of linear pol arization varies and orth ogonal mode jumps 
are present (|Karastergiou et 3111200911 . The bright pulses 
from PSR J0901—4624 show low levels of polarization, 
in contrast to its normal-mode emission, which at the 
same longitude shows a high level of linear polarization. 


There are two plausible interpretations for the different 
emission modes. The bright pulses could arise from a 
different plasma state or from a different region in the 
magnetosphere. 

Because the bright pulses originate close to the geomet¬ 
ric center of the pulse profile, they likely originate from 
low altitude, and is therefore unlikely to be depolarised 
by caustic effects. Therefore the low level of polarisation 
suggests that the emission arises from a different plasma 
state, for example when the polar gap is plasma-starved. 
The second state would have to produce normal-mode 
emission which is observed simultaneously to the bright 
pulses. 


7. CONCLUSIONS 

We have detected bright, narrow pulses from 
PSR J090I—4624. Similar to giant-micropulse emission 
observed in 4 other young energetic pulsars, the pulse 
energies follow a power-law distribution. Unlike other 
giant-micropulse emitters, the pulses originate from the 
center of the profile. Unlike the normal-mode emission 
for the pulsar, the bright pulses show low levels of to¬ 
tal polarization. These bright pulses could possibly be 
produced either when the magnetosphere is in a different 
state or from a different location in the magnetosphere. 
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